Parallel pumping of magnetic vortex gyrations in spin-torque nano-oscillators 
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We experimentally demonstrate that large magnetic vortex oscillations can be parametrically 
excited in a magnetic tunnel junction by the injection of radio-frequency (r/) currents at twice the 
natural frequency of the gyrotropic vortex core motion. The mechanism of excitation is based on 
the parallel pumping of vortex motion by the rf orthoradial field generated by the injected current. 
Theoretical analysis shows that experimental results can be interpreted as the manifestation of 
parametric amplification when rf current is small, and of parametric instability when rf current is 
above a certain threshold. By taking into account the energy nonlinearities, we succeed to describe 
the amplitude saturation of vortex oscillations as well as the coexistence of stable regimes. 
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Parametric excitations of magnetization oscillations 
have been extensively studied in the area of ferromag- 
netic resonance [T] in connection with the processes of 
spin- wave instability [2] and parallel pumping [5^. These 
phenomena are driven by the modulation of the oscilla- 
tion frequency of elementary excitations (spin-waves) as- 
sociated to the spatially uniform ground state obtained 
by a strong bias field. In the case of spin- wave instability, 
the spatially uniform mode is driven to large precession 
angle by radio-frequency (rf ) fields applied transversal to 
the bias field. The modulation of spin-wave spectrum is 
due to the nonlinear coupling of the uniform precession 
with spin-waves. On the other hand, in the case of paral- 
lel pumping, the rf field is applied along the direction of 
the bias field, and thus directly modulates the frequency 
of spin- wave modes j3j. 

In this letter, we investigate both experimentally and 
theoretically the implementation of the principle of par- 
allel pumping in spin-torque devices having a vortex in 
the free magnetic layer. It is an unconventional case of 
parallel pumping because the field used to modulate the 
frequency of magnetization oscillations is the orthoradial 
magnetic field generated by the injected current (usually 
referred as Oersted field). The frequency of this modu- 
lation is chosen to be twice the frequency of the lowest 
frequency mode present in a magnetic thin disk with vor- 
tex ground state, which is the translational motion of the 
vortex core [4]. 

The present investigation is relevant to the develop- 
ment and control of spintronics rf-oscillators in large ex- 
citation regimes [S] . These oscillators promise to have an 
important role in near future microwave communication 
technologies |6l[7] . A striking result in this areas has been 
the development of spin-torque nano-oscillators (STNO), 



tunable over a wide frequency range by the injected cur- 
rents both for uniform [8-10, and non-uniform [TIl413| 
magnetic configurations. The performances of these de- 
vices in terms of power and linewidth have been con- 
stantly improved over the last few years ^14-r,20j . In par- 
ticular, it has been recently demonstrated [2TI - [23] that 
Magnetic Tunnel Junctions (MTJ) pillars with vortex 
ground state [21H1S] allow to obtain signal with large 
output power (~ 1/iW) and very good coherence (~ 1 
MHz). It is crucial to gain understanding of the dy- 
namics of these oscillators in regimes where the vortex is 
driven very far from the equilibrium position as it is the 
case of the present study. 

Parametric resonance by means of parallel pumping 
in a spin-valve nanopillar has been previously studied 
by Urazhdin et al. (2!1\ in the more traditional setting 
where the magnetic free layer is uniformly magnetized 
and the pumping rf field, generated by an antenna, is 
applied along the spatially uniform bias field. In addi- 
tion, the somehow related phenomena of phase-locking 
of spin-torque auto-oscillations by an external source at 
double frequency, has been investigated both in the case 
of uniformly magnetized free layer |27| . and in the case 
of vortex state free layer [25]. 

Here, we investigate the parametric excitation in a 
vortex-MTJ subject to an external rf current with a 
frequency close to 2/o, where /o is the gyrotropic fre- 
quency of the vortex core with a dc-current below the 
threshold for self-oscillations (sub-critical case) used to 
reduce dissipative forces on the vortex. The device is 
subject to a strong out-of-plane static field, used to tilt 
the magnetization in the polarizing layer, which leads to 
an out-of-plane component of the vortex curling in the 
free layer 21- The parallel pumping is realized by mod- 
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FIG. 1: a-c) Colored maps of emitted power for the frequency 
of the gyrotropic core motion versus the external rf frequency 
measured (a) for a small (—18 dBm) and (c) a large rf current 
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spectral density versus frequency measured with an external 
rf current at frj = 275 MHz. 



ulating the gyrotropic frequency of the vortex through 
the rf variations of the Oersted field. The resulting gy- 
rotropic motion exhibits two specific features of paramet- 
ric resonance of weakly dissipative oscillators. The first 
effect, usually referred as parametric amplification (see 
e.g. [IS]), consists in the amplification of signals which 
contains the frequencies close to the natural oscillation 
frequency of the vortex. In our case, it manifests itself in 
the amplification of the thermally activated vortex mo- 
tion around its equilibrium position. The second effect 
is the parametric instability that takes place when the 
power of the rf excitation exceeds a certain threshold 
at which both the amplitude and the coherence of the 
oscillations strongly increase. Eventually saturation is 
reached due to nonlinearities. 

We perform our measurements on several circular 
MTJs of 500 nm diameter. The complete structure (with 
thickness in nm) is: PtMn(15) / CoFe(2.5) / Ru(0.85) 
/ CoFeB(3) / MgO(1.075) / NiFe(5) / Ta(7) / Ru(6) / 
Cr(5) / Au(200). The 5 nm-thick NiPe free layer presents 
a vortex magnetization at remanence. The MTJ resis- 
tance at the saturated state is i? ~ 45 with aver- 
age magnetoresistance Ai? ~ 8.5 at room tempera- 
ture. For the whole presented measurements, a dc cur- 
rent Idc = +3 mA and an out-of-plane field Hext,z = 4.48 
kG, are applied. In our convention, a positive current 
means electrons flowing from the free layer to the SAF 
stack. We emphasize that the injected dc current Idc is al- 
ways below the threshold Ith necessary to excite large am- 
plitude vortex-core sustained oscillations. Consequently, 
when no 7^/ current is applied, the small Idc induces only 
small oscillations of the vortex core and a very weak sig- 
nal is measured at /o — 138 MHz with maximum power 



about 3 nW/GHz/mA^ and a large linewidth (about 12 
MHz). These figures are typical of thermal- induced vor- 
tex core oscillations [21] |30l |3T] . The main features of 
the microwave signal are strongly modified when an ad- 
ditional current Irf is injected at a frequency close to 
2/o. In Fig. [T^, we first present the color-scale map of 
the power spectral density (PSD) as a function of the fre- 
quency of the external rf current at a small power equal 
to —18 dBm (corresponding to I^j — 0.8 mA). A signal 
emerges from the background level when 264 < f^f < 278 
MHz, which corresponds to a frequency around twice the 
natural frequency /o of the system. As shown in Fig[l]3, 
the spectrum measured for = 275 MHz can be fit- 
ted by a Lorentzian function and a maximum power of 
37 nW/GHz/mA^, about one order of magnitude larger 
than the thermal one, is obtained. The corresponding 
linewidth is also reduced of one order of magnitude to 
2.1 MHz. To our knowledge, these measurements are the 
first observation of parametric amplification for the vor- 
tex dynamics in the regime of sub-critical Idc currents 
and are consistent with the ones presented by Urazhdin 
et al. |27| for the case of uniform magnetization under a 
parametric rf field. 

The situation changes drastically when Irf > —18 
dBm. We enter in the parametric instability regime and 
a tremendous improvement for both the coherence and 
the emitted power is obtained. A striking example of 
such behavior is displayed in Figjl]: for large rf power 
equal to —12 dBm (corresponding to Irf = 1.59 mA). 
We find a large window of frf, i.e., between 265 and 
285 MHz, in which the emitted power associated to the 
vortex dynamics comes out from the background level 
and the spectral linewidth reduces significantly. In this 
range, the vortex frequency increases linearly with frf 
and is strictly equal to frf/'2. To demonstrate the dras- 
tic improvement of rf features in this specific regime, we 
plot in Fig. [Iji, the peak detected for frf = 275 MHz. 
The measured linewidth A/ of 49 kHz is indeed limited 
by the resolution bandwidth (RBW) used for such large 
range of measured frequency. Additional measurements 
with optimized RBW, allow us to extract a bottom limit 
value of the linewidth equal to 9 kHz. The maximum 
power reaches 23.6/iW/GHz/mA^ that is three orders of 
magnitude larger than the case of low rf power. 

In order to describe the phenomena observed in exper- 
iments, we consider the Thiele equation for the descrip- 
tion of translational motion of the vortex core [3TJ [5S] in 
the free layer of the MTJ device. This equation can be 
written as: 
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where X = {Xi,X2) is the vortex core position. The 
first term at the left-hand-side (LHS) of Eq. ([T]) is the 
gyrotropic term with G = —Bz^ttLMs/j 24J (7 the gy- 
romagnetic ratio, L is layer thickness, and Ms the satura- 
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tion magnetization). The unit vector is perpendicular 
to the plane of the layer and directed from the polarizer 
to the free layer. The second term at the LHS of Eq. ([I]) is 
the damping D = a'q\G\ where a is the Landau-Lifshitz 
damping, and 77 « 1 is a geometric factor. The first term 
at the right-hand-side (RHS) of Eq. ([T]) is the force asso- 
ciated to the gradient of the vortex energy W{X). This 
terms is given by the sum of the vortex magnetostatic 
energy Wmsi and the vortex Zeeman energy Woe associ- 
ated to the Oersted field which is proportional to the in- 
jected current density j = I/{ttR^) = >/ cos{2TTfrft)+jdc 
{R is the pillar radius). The expressions of these two 
energy terms, which introduce the main nonlinearity in 
the system along with the coupling with the RE exci- 
tations, are Woe — KejX'^ + KeJ^^/R^ and Wms — 
-(K™,i?V2) log(l-[X/2i?]2), where Aoe = OMfioMsRL, 
= -O.SAoe and k„,, = {10/9)^iqM^L'^/R [21 [21 US]. 
The term Fst,z — CzPzj x X is the spin transfer 
force responsible to compensate the damping forces. It 
depends on the perpendicular component of the polar- 
izer magnetization pz = p ■ Bz and on Cz — \G\'J(j/2, 
a constant measuring the spin-torque efficiency, where 
a — hP/{2\e\LMs), P is the spin polarization, h the 
Planck constant and e the electron charge. The term 
Feff at the RHS of Eq. ([T]) takes into account all addi- 
tional forces acting on the vortex. This terms includes 
spin-torque forces due to the in-plane component of the 
polarizer magnetization |23) . static in-plane forces act- 
ing on the vortex core (for example, the remaining stray 
field from the SAE polarizer) and extra forces due, for 
example, to defects [351 [33] . Such force Feff is assumed 
to be approximately independent of X. The main ef- 
fect of Feff under dc excitation is to produce a break of 
the rotational symmetry of the system which leads to a 
static shift of the core position from the MTJ center to 
the position Xq. The effect of this force is equivalent 
to an effective in-plane field Heff^xy Even if the precise 
determination of Heff^xy is not easy to obtain, from resis- 
tance versus in-plane field measurements, we estimate it 
to be few tens of Gauss, thus enough to break the system 
symmetry even at low injected rf currents. Einally, the 
strong out-of-plane external field is taken into account by 
an appropriate rescaling of the parameters entering in all 
terms of Eq. (|l]) I4,i23 . 

In order to study vortex oscillation around the dis- 
placed position, it is convenient to use the normalized 
core displacement Sx — {X — Xq)/R. After appropriate 
algebraic manipulations |34| . and by neglecting inessen- 
tial small terms, Eq. ([!]) can be rewritten as follows: 

^5x^Aq (t) -Sx + N (t, 5x) +v(t), (2) 
at 

where 

A (.^ _ f-di:!l{t) + CzPzJdc -^2{t) \ 

^l{t) -dW2{t) + CzPzJdcJ ' ^ ' 



d = D/\G\ is the normalized damping, Cz = Cz/\G\ is 
the normalized spin-torque efficiency and the two angular 
frequencies 

^iM^) = ^1,2 [1 + 91, 2>/ cos(27r/r/ t)] (4) 

are related to the quadratic approximation W{Sx) « 
[uji{t)Sxl + il!2{t)Sx2]/2 of the vortex energy around Xq. 

In Eq. UJi = n^e + ^'de^l, OJ2 = fide + ^U'^c^l, 

where a;o = \Xq\/R, fide = [nms + KeJ±)/\G\ and 
fl'^^ = {0.25Kjns + KeJct)/\G\. The frequency /o of vortex 
gyrations around Xq is given by the formula ujq = 27r/o = 
^LO\i>j2- The parameters q\ and (72 are related to the rj 
component of Woe'- 91 — (^oe + 3AQgXQ)/(a;i|G|) and 
92 = (\oe'rKeA)l{^^\G\) [3_4, . The term N (t,8x) ^ 
Nq {6x) cos{u!rft) in Eq. ^ contains all nonlinear terms 
in Sx. Finally, the last term v{t) = t>o cos(cJr/ i) is a 
forcing term independent of Sx, which comes from the rf 
components of Feff. 

The fact that in general uji (t) 7^ uj2 (t) is a consequence 
of the break of the symmetry due to the force Feff and 
it is indeed this break that controls the coupling between 
vortex oscillation and parametric excitations. The situ- 
ation is similar to the case of usual parallel pumping in 
rotationally symmetric systems. On that case, the cou- 
pling between the 2/o excitations and the magnetization 
oscillations is controlled by the ellipticity of the ampli- 
tude which is due to the dipolar fields generated by the 
spin- waves |3]. Instead, in our case parametric pumping 
of vortex dynamics takes place only if qi 7^ q2, which 
occurs when the vortex oscillates around a displaced po- 
sition. 

To treat analytically Eq. ([2]), it is important to no- 
tice that conservative terms in the RHS are expected to 
be the dominant terms in vortex core oscillations. The 
other terms act as perturbations on the system and thus 
induce only slow changes of the vortex oscillation. As 
a consequence, in order to describe this slow variation 
when ujrf ~ 2ujq, it is convenient to make the following 
(Van Der Pol type) change of variable Sx i~7> a: 

<5a:(t)= ('"°.'(:^'-/V?o\ ''^^^''\''''^'f:^'^^A-a{t) (5) 

^ ' yusm(uJrf t/2) COS{uJrf t/2) J w V / 

where u — {uji/uj2Y^'^ and cj^/ ~ "^"^frf- By substituting 
Eq. ^ into Eq. ^ and by taking the time-average of the 
resulting equation over one period of the rf oscillation, 
one arrives to the following equation 

da ^ /-duJo{a) + CzPzJde Suioia) - Jrf'^ \ ■ a (6) 
dt \ Suio{a)~jrfT ~duio{a) + CzPzJdc) 

where SQo^a) ~ iOQ{a) — uJrf/2 is the nonlinear detuning 
parameter and 

uJo{a) =ujQ + VoeidcO^ + rirns[l/(l - (a/2)^)] , (7) 

with flms = >^ms/\G\ and Voe — ^oe/\G\. The parame- 
ter which control the coupling with the rf excitations is 
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FIG. 2: a) Phase diagram in the (Z^/, jrf )-plane with indica- 
tion of region of parametric instability, coexistence of stable 
regimes and parametric amplification, b) Analytically com- 
puted amplitude response versus frf, according to Eq. (Ill 
with Irf — 1.78 mA which corresponds to —11 dBm rf power. 
The values of the parameters used are: P = 0.225, Idc = 2.5 
mA, a — 0.01 (LL damping) and xq = 0.3. 



FIG. 3: a) Vortex frequency / versus r/-frequency frf- b) 
Power density versus r/-frequency frf. c-d-e) Power density 
versus frequency / for three values of frf : 275, 295 and 305 
MHz. Note that these measurements have been made with an 

4.48 kG, a subcritical dc current 



out-of-plane field Hperp = 
Idc = 2.5 mA and a rf current /, 



rf 



-11 dBm. 



jrfT = Wo (92 — 9i)/4 which is zero if there is no shift of 
the vortex, i.e. qi =52- 

Note that, in the condition under study, i.e. Idc < hh, 
there is no sustained oscillation and it is the rf time- 
varying Oersted field contribution which is responsible to 
destabilize the vortex core from its equihbrium position. 
The stability of small magnetization oscillations around 
the static equilibrium position can be studied by consid- 
ering the linearized version of Eq. ^ around a = 0. The 
eigenvalues of such linearized equation are given by 

Ai,2 = -dujQ + c^pzjdc ± \/jlf^^ - ('^o - '^rfl^Y ■ (8) 

The stability is controlled by the sign of the real part of 
Ai^2 and the instability boundary in the (/r/, >/) -plane 
correspond to the condition of vanishing of one eigen- 
value: 



2 p2 



(9) 



This is the equation of the hyperbola centered at the 
point (wr/ = 2ujQ,jrf = 0) and with asymptotes along 
the lines >/ = ±(a;r//2 — wo)/r (see Fig. [2^). From 
Eq. ^, one derives that the threshold value of jrf such 
that instability occurs is 



th 

Jrf 



(10) 



For the particular parameter choice of Fig. [2^, j*y — 1 
mA. Below j'j we find the very same behavior observed 
in Fig. [T^-b. In general, outside the hyperbola the real 
parts of Ai 2 are negative and small vortex oscillation are 
stable. Nevertheless, if one consider linear response of the 
system to noise, one may observe two different qualitative 



behaviors. At relatively large detuning, the eingenvalues 
are complex conjugates and the frequency response of a 
has peaks at the frequencies /r//2 — /o and /o — /r//2. 
By using Eq. (§, i.e., going back to the stationary frame, 
one can infer that thermally driven vortex dynamics have 
power spectral density doubly peaked at /q and frf — fo- 
The behavior changes close to 2/o (blue region in Fig. [2^) 
when the eigenvalues become both real. The frequency 
response of a is centered at zero, and thus, by using again 
Eq. (§, one can conclude that the power spectral density 
of thermally driven oscillations is centered at /r//2. The 
amplitude of the linear response increases progressively 
as {frf , Jrf) approaches the instability boundary and the 
signal linewidth shrinks to zero. This behavior is the 
so-called parametric amplification. In our experimental 
investigations, the evidences of parametric amplification 
are seen in the increase of thermally driven vortex oscil- 
lations which are precursors of parametric instability, see 
Fig§a-b). 

Inside the instability region (red region in Fig. [2|, the 
equilibrium position of the vortex becomes unstable and 
the study of the system response requires the inclusion of 
nonlinear terms. The nonlinear regime after instability 
can be found by searching for nonzero equilibria of the 
general equation Eq. (|6| which, by taking into account 
Eq. ([5]), correspond to steady state of the vortex at a 
frequency /,//2. This solution can be found by imposing 
that the determinant of the matrix at the RHS of Eq. ^ 
is zero, which leads to the equation: 

jlfV"^ - [wo(a) -ujrf/2f = [du)n{a) - c^PzJdcf ■ (H) 

The nonlinear response of the system can be then ob- 



tained by keeping j^/ fixed and interpreting Eq. (Ill as an 
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implicit relation between a and w^/ ■ The result of this an- 
alytical computation for the particular case of Irf = 1.78 
mA (dashed black line in Fig. [2^) is shown in Fig. [2J3. 
Solid lines indicate stable regimes, while the dashed lines 
indicate unstable regimes. In the white region only one 
state exists and it is stable (a = 0). In the red region, 
parametric excitation takes place and the amplitude of a 
start to increase while the state at a = becomes unsta- 
ble. In the region of coexistence (green region in Fig. [2]) 
the stable large amplitude regime coexists with the stable 
equilibrium position of the vortex. 

To verify the validity of the analytical model, notably 
the prediction of the coexistence (green region in Fig. [2]), 
we present in Fig. [3] a serie of measurements performed 
on an other junction from the same wafer, in which the 
vortex gyrotropic dynamics has a few MHz difference for 
the thermal resonant frequency. The measurements have 
been performed with Hperp = 4.48 kG, Idc = 2.5 mA and 
7^/ = -11 dBm. 

In Fig. |3^-b, we plot the vortex frequency / and the 
power of the emitted signal as a function of frf- At low 
frequencies {frf < 275 MHz), we detect only a small peak 
around /o corresponding to the thermally excited signal. 
At frf = 275 MHz, a second peak with a very small 
power (note the logarithmic scale in Fig. [3]d) centered at 
frf — fo appears. This situation with two small peaks, 
whose power increases with frf lasts until 285 MHz. A 
typical spectrum in this regime is shown in Fig.|3j:. From 
frf = 285 until 305 MHz, the spectra contains a main 
peak with large power and two side bands. The anal- 
ysis of these thermally excited sideband signals in the 
parametric regime is out of the scope of this paper and 
for sake of clarity, we have not reported theirs values 
in Fig. |3^-b. In this rf frequency range, the vortex fre- 
quency is locked to the source signal, i.e., frf/2, and thus 
evolves linearly (see Fig. [3^). As shown in Fig. ^p, the 
peak power increases strongly to reach a maximum of 
« lO^iW/GHz/mA^ at frf = 295 MHz. The linewidth 
obtained at this rf frequency is only a few tens of kHz. 
Then, between 295 and 300 MHz, the power spectral den- 
sity stagnates, meaning that the amplitude of the vortex 
oscillation saturates. Increasing further frf above 300 
MHz, the power density start to decline. Moreover at 
frf — 305 MHz, we detect again a second peak centered 
at fo (see Fig.[§). 

As expected by the analytical calculations (see 
Fig. [2]d), the two peaks, one at /r//2 and one at /q, 
detected in the region between 305 and 315 MHz, cor- 
respond to the two states predicted in the region of coex- 
istence. Notably, the main features of the power density 
evolution (see Fig. [sja) is in very good agreement with the 
evolution expected by the theory (see Fig. [2]d) : the oscil- 
lation amplitude grows monotonously from 100 to 10000 
nW/GHz/mA^ between 282 and 295 MHz. Conversely, 
we find that the opposite reduction of power occurs very 
sharply in few MHz (from 300 to 305 MHz). This behav- 



ior corresponds to the soft-hard regime discrimination al- 
ready observed for the case of uniform magnetization [2 7) . 
Note that thermal effects, that are not considered in the 
model, result in a blur of the transition between the dif- 
ferent regimes, and consequently allow some thermally 
induced transitions from one state to the other. Further- 
more, we believe that another impact of thermal energy 
is to avoid the experimental observation of the hysteresis 
effect predicted in Fig. [2]d. Complementary experiments 
at low temperature and analysis of the parametrically 
driven vortex dynamics in time-domain might be helpful 
to address more precisely the impact of temperature on 
the non-linear vortex dynamics. 

In summary, we have presented a comprehensive inves- 
tigation of parametric excitation of vortex dynamics in 
a MTJ based spin transfer oscillator. Moreover, we pro- 
pose an analytical model to predict the phase diagram of 
our vortex system in presence of an external rf current 
at about twice the natural vortex frequency. We report 
the first observation in spin torque devices of the differ- 
ent parametric regimes of amplification and instability. 
Finally, we demonstrate both experimentally and theo- 
retically, the coexistence of two parametric states of the 
vortex dynamics, evidencing the important role of non 
linearities. The parametric excitation in vortex based 
devices might be used for highly efficient rf detection 
or low noise amplification using the specific potential of 
these devices in tuning their rf features through both the 
rf and da current. 
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